o 
o 



o 



1 



On the application of differences in intrinsic 

fluctuations of Cherenkov light images for 

separation of showers. 



V. V. Bugayov^'^ A. V. Plyasheshnikov ^ V. V. Vassiliev^ 

T. C. Weekes^ 



^Department of Physics, Altai State University, Dimitrova 66, Barnaul, Russia, 

656099 

^Fred Lawrence Whipple Observatory, Harvard-Smithsonian CfA, P.O. Box 91, 
^ . Amado, AZ 85645-0097 

> 

Q> , Abstract 

in 



The sensitivity of ground-based imaging atmospheric Cherenkov 7-ray observatories 

depends critically on the primary particle identification methods which are used to 

retain photon-initiated events and suppress the spurious background produced by 

^S^ ■ cosmic rays. We suggest a new discrimination technique which utilizes differences in 

*^| the fluctuations of the light intensity in the images of showers initiated by photons 

O ' and those initiated by protons or heavier nuclei. The database of simulated events 

for the proposed VERITAS observatory has been used to evaluate the efficiency of 

j^ ! the new technique. Analysis has been performed for both a single VERITAS imaging 

telescope, and a system of these telescopes. We demonstrate that a discrimination 

efficiency of > 1.5 — 2.0 can be achieved in addition to traditional background 

rejection methods based on image shape parameters. 
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1 Introduction 



Gamma-ray astronomy in the energy range above 100 GeV has made dra- 
matic progress in the last decade due to rapid development of the atmospheric 
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Cherenkov technique at a number of ground-based 7-ray observatories. Detec- 
tion of Cherenkov radiation from atmospheric cascades initiated by high en- 
ergy particles has been considered as a potentially promising method for high 
energy 7-ray astronomy since the fifties [1] due to the large collecting area of 
the technique, ~ 5 x 10'^ m^. Only recently, however, has it become possible to 
overcome the main factor limiting the sensitivities of ground-based observa- 
tories which utilize this method. This is the high rate of the cosmic ray (CR) 
background events. Introduction of the two techniques, imaging and stere- 
oscopy, pioneered by the Whipple and HEGRA collaborations respectively 
[2,3], has made a drastic improvement in the suppression of the background, 
and has made possible the reliable detection and spectroscopy of several galac- 
tic and extragalactic sources [4-7] . Recent reviews of the development of high 
energy ground-based 7-ray astronomy can be found in [8,9]. 

Methods for the efficient discrimination of photon and hadron initiated events 
have been derived from the differences in the intrinsic properties of the Cherenkov 
radiation from purely electromagnetic and hadronic cascades. The most promi- 
nent of these is the substantially wider distribution of the arrival directions 
of photons in hadronic showers due to an additional angular scattering of the 
secondary products during hadronic interactions. This effect has been exten- 
sively studied in Monte-Carlo simulations of the cascade development in the 
atmosphere [18,17], and the rejection criteria for the cosmic ray background 
has been derived in the form of the so-called "shape cut" [10,11]. Increased an- 
gular resolution of ground-based 7-ray observatories together with this means 
of identification of the primary photon are the most important components of 
the recent discoveries and advances made in the highest energy 7-ray astron- 
omy. 

The confidence of discovery of a 7-ray source in the case of a background 
dominated data sample is defined as 



ij = %i5£, (1) 

where A^on is the number of detected events produced by a hypothetic source 
for a given exposure, and Nos is the number of background events for the 
same observation time. This quantity refiects an excess of events due to the 
presence of a source measured in units of background standard deviation shown 
in the denominator for a Poisson distribution. It is generally accepted that the 
existence of a source is established ii R> 5. 

By application of the background rejection procedure one increases the confi- 
dence of the source discovery due to elimination of the signal fiuctuations con- 
nected with the background. The efficiency of the discrimination technique 
used to identify primary photons and reject cosmic ray events is expressed 



then as 



V = ^. (2) 

where k^ and k^ are the probabihties of acceptance of photon and background 
events after the apphcation of selection criteria. Thus the rj factor is an en- 
hancement of R indicating that rj^ times less exposure is required for a source 
discovery when a new discrimination technique is applied. This is true in the 
case of background dominated regime, when N^^ — Nos < ^off- On the con- 
trary, in the statistics dominated regime (A^on — ^ofr >> ^ofr), background 
rejection procedures with k^ remarkably smaller than 1 should be avoided. 

In this paper we explore a new discrimination technique that is characterized 
by efficiency rj = 1.5 — 2 in addition to the already existing methods of back- 
ground suppression achieved by imposing shape and orientation cuts. This 
technique is based on the differences in intrinsic fluctuations of Cherenkov 
radiation produced in the pure electromagnetic and hadronic cascades. We 
adapt the maximum likelihood approach based on an analysis of the x^ func- 
tional to establish discrimination criteria to distinguish showers initiated by 
7-rays from those initiated by CR nuclei. This idea was induced originally 
by the study of the mass composition of primary cosmic radiation [12] where 
similar methods were proved to be very successful. 



2 The intrinsic fluctuations of the image 



The imaging atmospheric Cherenkov telescope (lACT) consists of a set of mir- 
rors on a single mount with a two-dimensional array of close-packed hexag- 
onal photomultiplier tubes (pixels) in the focal plane. Therefore information 
about an event is recorded in the form of a two-dimensional image describ- 
ing the distribution of the light intensity in the focal plane. Ultimately, the 
Cherenkov light image is a set of ADC counts corresponding to all pixels. 
But for simplicity we deal with the number of photoelectrons emitted from 
the pixels' photocathodes. Hence the Cherenkov light intensity in the focal 
plane is characterized by the continuous two dimensional distribution of pho- 
toelectron density, p{x,y). Each point {x,y) on the focal plane corresponds 
to a certain arrival direction of photons relative to the pointing of the tele- 
scope mount. The total number of photoelectrons (pes) in the Cherenkov light 
image, SIZE, is given by 
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In our further calculations we make frequent use of the normalized to unity 
distribution of photoelectrons 



f{x,y) = SIZE-'p{x,y), (3) 

which is equivalent to the probability density function (pdf). By translation 
and rotation of the reference frame it is always possible to choose a new 
coordinate system {x^y) in which the first moments of f{x,y) are equal to 
zero and the matrix of second moments is diagonalized. The diagonal elements 



fx^f{x,y) dxdy = LENGTH^ 
fy^f{x,y) dxdy = WIDTH\ 



which characterize the extent of the Cherenkov light distribution along the 
major and minor axes of the image, are often used as the image shape param- 
eters. 

Let us denote also the marginal distributions 



fL{x)= j f{x,y) dy, 
fw {y)= f (x, y) dx, 



where indices L and W correspond to the major and minor axes of the image 
respectively. We estimate /^ (x) and fy/ {y) by the grid functions 






'L,k = J II {x) dx, Xk+i = Xk + Ax; 

Xk 
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Nw,k = fw iy) dy, yk+i = yk + Ay; 
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defined at {xk, yk',k = l,..., K}. 

Suppose that we have two sets of shower images, initiated by photons and 
cosmic rays, for which the mean values, N^l , N^f], and variances, o"^ (^), a^ (p) 

t,k t,k 

of the grid functions are found for each k,t G {L, W}. For each image we 
construct a pair of functionals Xt 
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SO that the mean values of Xt fo^ 7-showers and background events are equal 
to: 



Xt 
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The Xt^^' consists of two terms. The first term reflects the differences in the 
fluctuations of the light intensity of images initiated by 7-rays and those ini- 
tiated by CR particles. The second one characterizes the differences between 
the average shapes of distributions of the photoelectron density in the images 
of the two data sets. 



Strickly speaking, two different effects influence the value of the "fluctuation" 

term of Xt ■ The first of them is connected with random variations of the 
angular size of the image described, for example, by WIDTH and LENGTH 
parameters. The second effect is the irregularity of the light distribution inside 
the image itself. We say in this case about the intrinsic fluctuations. Even in 
the case of 7-ray and proton induced images having the same angular size 
the first of them has a more smooth and regular structure and, therefore, 
a smaller value of intrinsic fluctuations. Large fluctuations in the number of 
secondary particles created during the hadron multi-particle production is the 
main reason of difference in intrinsic fluctuations of 7- and p-induced images. 



In this work we attempt to make use of the contribution to Xt from the 
first ("fluctuation") term because we expect that the second term (as well as 
the fluctuations connected with variations of the angular size of the image) is 
rather small for images which passed the image shape selection criteria. 

Keeping in mind the argumentation presented above, one can hope to use 
selection criteria of the form 



2 ~2 

Xl < Xl, 



2^-2 

Xw ^ Xw 



(6) 



to reject CR-induced events and retain genuine photon-initiated showers. For 
a single telescope, two constants Xl ^^^ Xw should be optimized to maximize 
the signal to noise ratio (see formula (1)). For a system of telescopes a set of 
similar criteria can be utilized 



xL 



<xl, 



2 ^ ~2 

Xw,i ^ Xwy 
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trig. 
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in which condition (6) is satisfied simultaneously for all triggering telescopes. 
We consider a slightly more relaxed background rejection criterion 



1 -^trig 1 -W^trig 

-L X — ^ o ^o r 



E2 ^ ~2 V^ 2 ^ ~2 



which seems to retain more photon-initiated events and produce a higher 
efficiency factor (2). In the case of a single telescope, both criteria (6) and (8) 
are identical. 

The success of a selection criterion depends dramatically on how the data sets 
of photons and CR nuclei have been selected. In order to increase the selec- 
tion efficiency one should consider the smaller regions of the primary photon 
parameter phase space and formulate selection criteria for each of them. We 
expect that the shower parameters which affect background rejection most 
strongly are the primary energy, E and the impact parameter, r (the shortest 
distance between the telescope and the cascade core). In order to make the 
X distribution narrower and thus improve the discrimination ability of the 
method, we will evaluate the dependence of the grid functions of photoelectron 
density on these parameters 



N,,, = N,,,iE,r), Nil^ = Ni]]iE,r), a^(,) = a^^ (i?, r) . (9) 

t,k t,k 



We would like to note that the x^-based method has been used in the past to 
deduce parameters of the primary particle [13,14]. In our work, however, we 
do not explicitly assume that the ffuctuations in the images of 7-rays are of a 
Poissonian nature. 



3 Data sets of shower images 



To investigate the discrimination efficiency of the method we have simulated 
two data bases of photon and proton initiated showers. This task has been 
accomplished by utilizing ALTAI computation code [15] with input data sum- 
marized in Table 1. The geometrical layout of the VERITAS array of telescopes 
is shown in Fig. 1. Although, both proton and photon induced cascades were 
simulated as vertical showers, images of cosmic rays were also sampled isotrop- 
ically on the focal planes of the telescopes near the vertical direction. 44,000 
and 92,000 showers were simulated for photons and protons respectively. The 
position of the shower core was distributed uniformly inside a circle with ra- 
dius of 300 m centered on the central telescope of the VERITAS array. A 






Fig. 1. Layout of VERITAS telescopes. 

power law differential energy spectra of primary particles was assumed with 
exponent —2.6 for photons and —2.75 for protons. Photons were sampled 
from the energy interval 0.05 — 10 TeV and protons were distributed between 
0.1 -20TeV. 

The contribution of the night sky background was taken to be one photoelec- 
tron per pixel. We applied a rather severe cleaning algorithm which insured 
that only bright images with light intensity well above this background were 
included in our analysis. Shower images were processed in two stages. First, 
all pixels with less than five photoelectrons were excluded from further consid- 
eration. Then, only pixels which contained more than nine photoelectrons or 
which had a neighbor pixel with a nonzero content were assumed to be image 
pixels. In addition, only images with SIZE greater than 100 photoelectrons 
were accepted for the analysis. The telescope was considered to be triggered 
when at least two pixels from an internal zone of the multichannel camera 
(331 inner pixels) detected more than 20 photoelectrons. Finally, to ensure 
the high quality of images, only telescopes with maximal light intensity of the 
image in the internal zone were accepted. 



Table 1 

Basic parameters of VERITAS telescopes 

Number of telescopes 7 

Altitude above the sea level, [m] 1500 

Area of the mirror, [m^] 78 

Number of pixels 547 

Pixel size, [degree] 0.15 

Field of view, [degree] ~4.0 



4 Determination of the impact parameter 



The efficiency factor, rj, of the discrimination teclinique increases if the de- 
pendence of the image photoelectron density on the energy and the impact 
parameter of a shower is accounted for (formula(9)). Therefore, it is necessary 
to have some method of determining these quantities. In this work we have 
investigated three possibihties for estimating impact parameter of the primary 
particle: 

• VERITAS operates as a system of telescopes. The impact parameter is de- 
termined by a simple geometrical event reconstruction method (see, e.g., [11]). 

• VERITAS operates in a single telescope mode. 

• The impact parameter is determined on the basis of the DIST (distance) 

parameter of the image. 
■ The impact parameter is determined by means of the ELLIPT (elliptic- 

ity) parameter of the image. 

Stereoscopic observation of an event by at least two telescopes allows unam- 
biguous geometrical reconstruction of the shower core location in space. If the 
number of telescopes which detected an event is larger than two it is possible 
to increase the accuracy in determining this parameter. In this work we use 
as an estimate of the shower core location, a weighted average of the core 
coordinates derived from the data of each pair of triggered telescopes [13]. 

For a single telescope, determination of the impact parameter utilizes its cor- 
relation with such characteristics of the image, as ELLIPT or DIST [16] 
defined as 

ELLIPT = ^^^^^^ ~1, DIST = ^xl + yl (10) 



where Xc and yc are coordinates of the image centre of gravity in the reference 
frame with an origin at the source location. Unlike the DIST parameter, 
ELLIPT provides an opportunity for determining the impact parameter not 
only for a point-like, but also for an extended 7-ray sources when the arrival 
direction of a photon is not known. To derive our estimate of impact parameter 
in the case of a single telescope observations, we tabulated mean values of 
ELLIPT and DIST for a set of values of impact parameter and image SIZE. 
The required estimate of r was then calculated by backward interpolation. In 
fig. 2 we show the dependence of mean DIST on the impact parameter for 
several intervals of SIZE. Larger values of SIZE correspond to curves with 
larger values of DIST. 

In fig. 3 and tab. 2 we show the accuracy of determination of the impact pa- 
rameter by three different methods. 
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Fig. 2. Dependence of the DIST parameter on the impact parameter for 7-showers. 
The curves are calculated for SIZE intervals with the centres at 220 (bottom curve), 
840, 3200 pes. (top curve) respectively. 
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Fig. 3. The probability distribution for the error of determining the 7-shower impact 
parameter. The SIZE of the images is larger than 150 pes. 

The accuracy of determining the impact parameter depends strongly on the 
lower bound of the SIZE of images included in the analysis. The smaller 
the SIZE of the shower image, the stronger the effects of fluctuations in 
image photoelectron density, and the stronger the influence of the night sky 
background. A large difference between the errors of impact parameter de- 
termination for photons and protons in the case of a system of telescopes is 
due to the inherent differences in the electron density distribution between 
the images of showers initiated by these primaries. While the distribution of 
Cherenkov light is smooth and compact for 7-rays, for protons it is fuzzy and 



more extended. As a result the proton induced images (and their basic param- 
eters determining the reconstructed value of the shower impact parameter) are 
more significantly distorted in the process of the cleaning procedure aimed at 
the elimination of the night sky background and consisting in an exclusion 
from consideration of pixels with low magnitudes. 

In the case of the photon primary the method of impact parameter determina- 
tion based on the DIST parameter is substantially more accurate than the one 
where the ellipticity is in use. Perhaps, this can be explained by the larger fiuc- 
tuations of the ELLIPT parameter, or in other words by weaker correlation 
of ELLIPT with shower core location. The "ELLIPT" based method indi- 
cates a substantially smaller error on the impact parameter determination for 
7-induced air showers relative to proton initiated ones, which can also be ex- 
plained by the larger fluctuations of the light intensity in CR-induced images. 
For isotropically distributed cosmic rays one should not expect any correlation 
between the impact parameter of the showers and the DIST value. The error 
on the determination of the shower core location of 53 m, found in this case, is 
finite because of the limited telescope field of view and the high directionality 
of the distribution of Cherenkov photons from atmospheric cascades. 



5 Event selection based on the shape parameters 



To estimate the discrimination efficiency of this new technique, we applied a 
preliminary event selection on the basis of the shape parameters {WIDTH and 
LENGTH) thereby excluding any possible correlation between the traditional 
and proposed event classification methods. In this work we adapted a selection 
criterion in the form 



WIDTH 

W('r){SIZE,^ 



<W,,, ^=l,...,iVtrig (11) 



where W^'^^ {SIZE, r) is the mean value of WIDTH for 7-showers detected 
with a given image SIZE and estimated impact parameter, r, (fig. 4). Ap- 

Table 2 

The value of the impact parameter determination error [m]. 50% of showers have 
an error smaller than the value presented in the table. The SIZE of the images is 
larger than 150 pes. 

7 P 



System of telescopes 5 32 

Single telescope {DIST) 11 53 

Single telescope {ELLIPT) 27 68 



10 



plication of this criterion, commonly known as the scaled width method [11], 
turns out to be particularly effective as it accounts for the changes of image 
WIDTH caused by the differences in primary photon energy and the possible 
variations in the shower core location. We excluded LENGTH as a discrimi- 
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Fig. 4. W'^ = f{SIZE,r). The curves are calculated for SIZE intervals with the 
centres at 220 (bottom curve), 430, 840, 1600, 3200, 6200 pes (top curve), respec- 
tively 

nation parameter while using Wsc] two-dimensional discrimination on the basis 
of LENGTH and WIDTH yields only negligible increase in r] [11]. 




Fig. 5. Dependence of background discrimination efficiency rj, true classification 
probability ( k^) of 7-rays, and false classification probability {np) of protons as 
functions of the Wgc cut. In this example, VERITAS operates as an array of tele- 
scopes. 

In fig. 5 the dependence of background discrimination efficiency on Wsc cut 



11 



(11) is shown. The dependence of t] on the true classification probabihty of 
7-rays is shown, for different event reconstruction methods, in fig. 6. 
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Fig. 6. Dependence of the background discrimination efficiency ij on the true clas- 
sification probability k^ of photons. 



6 Optimization of the discrimination technique 



Fig. 7 depicts two histograms: the ratio of the mean photoelectron density 
along the major axis {t = L) for images initiated by photons and protons 
(1), and the ratio of their standard deviations (2). The SIZE parameter was 
limited within interval from 100 to 200 photoelectrons, and the actual shower 
impact parameter was restricted to the 50 — 75 m interval. The values of these 
intervals were choosen as small as the statistics allows in order to reduce the 
fluctuations connected with random variations of the primary energy and of 
the impact parameter. 

It can be seen that for the given conditions, the contribution to the x^ func- 
tional (Eq. 5) from the first, "fluctuation", term exceeds the second, "mean 
difference" , term by almost a factor of ten. This confirms our assumption about 
the potential of the Cherenkov light intensity fluctuations for discriminating 
proton and 7-induced images. 

The bin sizes. Ax = 0.2° and Ay = 0.06°, used to construct the grid photoelec- 
tron density functions, were optimized to achieve the highest discrimination 
efficiencies. They are chosen as a balance between the necessity to resolve fea- 
tures in the shower images on the scales smaller than LENGTH and WIDTH 
and, at the same time, not to make fluctuations in the bins dominated by Pois- 
son statistics. In addition, a special analysis has shown that the x^-technique 
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Fig. 7. NLkl-^Lk i^) ^^^ ^ n(i'> / ^ N^p'i (^) ^•^'^ events surviving the scaled width cut 

{Wsc = 1-4, DIS'T-approach) . A single VERITAS telescope is considered. Bin size 
Ax = 0.2°. 

exhibits no essential sensitivity to the bin size if the value of this size does not 
exceed considerably the pixel size (0.15") of the VERITAS telescopes. 

The comparison of the probability distributions xi and Xh' fo^ 7 and p-induced 
events, for a single telescope, is shown in fig. 8. Tab. 3 summarizes the photon 
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Fig. 8. Probability distributions for Xl (*-"^ ^^^ ^^^^) ^^"^ Xw (^^ ^^^ right) for a 
single VERITAS telescope after Wsc = 1-3 cut is applied. Determination of the 
impact parameter is based on DIST parameter. 

selection efficiency after application of the scaled width cut, {k^WiVw), and 
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after discrimination of background using the Xl ^^^ Xw parameters (eqn. 8), 
( K-y^, r]y.). These calculations were performed for the full VERITAS array. The 
first column indicates the value of VK,c cut. The total discrimination effect is 
denoted as («;^toti ''7tot)- '^7x corresponds to the value of xi ^-nd Xw ^^^^ ^^ 
which discrimination efficiency factor, r], is maximal provided that k^^ > 0.7. 

Table 3 
Summary of the background rejection efficiency for VERITAS. 

Wsc K^W VW K^x Vx '^7tot ??tot 

1.15 0.52 4.11 0.81 1.41 0.42 5.80 

1.20 0.60 4.08 0.74 1.43 0.44 5.83 

1.25 0.67 3.90 0.79 1.58 0.53 6.16 

1.3 0.72 3.87 0.82 1.66 0.59 6.42 

1.4 0.79 3.55 0.70 1.84 0.55 6.53 



It is evident from the table that the largest value of the total efficiency fac- 
tor, r^tot; is achieved for Wgc larger than the value which provides maximum 
discrimination power of photon selection based solely on the scaled width cut. 
This is also correct for event reconstruction using a single telescope. Selection 
criterion utilizing only the Xl cut demonstrates a discrimination efficiency fac- 
tor ~ 80 — 90% of the one obtained with the use of both x^ cuts. An increase of 
Wsc provides a moderate growth of photon discrimination efficiency when Xw 
is sequentially applied, while the effect of a xi selection depends very weakly 
on the scaled width cut. 

The plots of discrimination efficiency of x^ technique as a function of true 
photon classification probability, k^, are shown in fig. 9 for various VERITAS 
event reconstruction methods. The scaled width cut was optimized to provide 
maximal total discrimination efficiency when k,^ is fixed. Curves (1) and (2) in 
this figure compare the factor r] for a x^ cut in the form given in (8) and (6) 
when the latter is applied to each VERITAS telescope in the event trigger. It 
can be seen that the discrimination efficiency of photon selection is higher for 
the case which uses the mean value of x^- 

In tab. 4 we present the true photon classification probability and discrimina- 
tion efficiency of the x^ technique when various lower bounds on image SIZE 
are included in the analysis. The dependence of these quantities on the lower 
limit of the estimated energy of the primary particle is shown in the tab. 5. 
For both tables we tuned the scaled width cut to optimize the total efficiency 
factor. Primary energy was estimated by backward interpolation with respect 
to image SIZE and estimated impact parameter, r, using simulated data base 
of photon induced showers. 
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Fig. 9. Dependence of x^ selection efficiency rj (k^) for different event reconstruction 
methods. For array of telescopes, curve (a) corresponds to x^ selection criterion (8), 
and curve (b) is obtained for criterion (6) applied to each telescope participating in 
the VERITAS trigger. 

Table 4 

Discrimination efficiency of x^ photon selection technique ( W^c cut is tuned to 
maximize ??tot)- Different image SIZE lower bound for analyzed events is shown. 



System 



Single(D/5T) 



SingleiELLIPT) 



S 100 150 200 300 100 150 200 300 100 150 200 300 
0.75 0.70 0.71 0.71 0.73 0.70 0.73 0.70 0.79 0.71 0.70 0.70 
2 1.50 1.84 1.94 2.14 1.43 1.81 2.34 4.11 1.22 1.37 1.56 1.95 






^x 



Table 5 

Discrimination efficiency of x^ photon selection technique as a function of lower 
bound on event estimated energy to be included into the analysis. 



System 



Single (DIST) 



Single (ELLIPT) 



E^cst, [TeV] 0.2 0.3 0.5 1 0.2 0.3 0.5 1 0.2 0.3 0.5 1 






0.72 0.73 0.70 0.71 0.70 0.70 0.71 0.70 0.72 0.71 0.71 0.70 
1.91 2.00 2.19 3.33 1.62 1.97 2.31 3.92 1.52 1.65 1.76 2.47 



7 Conclusion 



In this work we have examined the feasibility of improving photon selection us- 
ing differences in the fluctuations of the Cherenkov light distribution between 
the images of 7 and p -initiated showers. Three different event reconstruction 
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methods which are hkely be utihzed in the next generation of ground-based 
7-ray observatories, such as VERITAS, have been studied. Two of them are 
apphcable to observations with a single telescope and one uses the stereoscopic 
capability of VERITAS array. 

Optimization of our technique indicates that the proposed method provides an 
additional increase of the total discrimination efficiency by a factor > 1.5 — 2.0 
retaining the survival probability of 7-ray initiated events at the level larger 
than 0.7. That is true for both the telescope array and the single telescope 
(with estimation of impact parameter by the DIST-hased method). The 
achievable value of the total discrimination efficiency is equal to 6.5 and 
3.6 respectively. When a single telescope event is reconstructed utilizing the 
ELLIPT parameter, our approach provides an additional discrimination effi- 
ciency factor > 1.4. This is particularly important for a detection of extended 
sources or sources whose position has not been accurately identified . The 
discrimination efficiency of the proposed method increases rapidly with the 
energy of the primary photon. 
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